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1. Introduction 

We have developed a bulk PBL model with a simple internal vertical structure and a simple 
second-order closure, designed for use as a PBL parameterization in a large-scale model. The 
model allows the mean fields to vary with height within the PBL, and so must address the vertical 
profiles of the turbulent fluxes, going beyond the usual "mixed-layer" assumption that the fluxes of 
conservative variables are linear with height. This is accomplished using the same convective mass 
flux approach that has also been used in cumulus parameterizations. The purpose of this brief 
paper is to show that such a mass flux model can include, in a single framework, the 
"compensating subsidence" concept, downgradient mixing, and well-mixed layers. 


2. Generalized mass flux model 


The "convective mass flux" concept introduced by Arakawa has been used in a variety of 
modeling and observational studies, most directed either at deep cumulus convection or at 
stratocumulus convection; the connections between the two applications have not been made 
previously. 

For an ascending region with fracitonal area o, the conservation equations for mass and an 
arbitrary scalar \\f can be written as 
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for the ascending parcels. Similar equations can be written for the descending regions. Here the 
generic variable \\f represents an intensive thermodynamic variable such as the potential temperature 
or the the mixing ratio of water. Other notation is as follows: p is the density; V is the horizontal 
velocity vector; co = -pgw, where w is the vertical velocity, and g is the acceleration of gravity; and 
p is pressure. The air entering rising parcels has been assumed to have the average properties of 
the sinking parcels, and vice versa. Area averages satisfy 

C~~)=( W( 

(2.3) 

where an overbar denotes an area average (over a grid box, say); this notation will be used only 
where necessary to avoid confusion. Subscripts u and d denote upward and downward moving 
parcels, respectively. We can interpret p-‘ and v> as time scales for mass flow from downdrafts 
into updrafts, and from updrafts into downdrafts, respectively. The vertical turbulent flux of \|/ 
due to small eddies is denoted by f ¥ . The small-eddy fluxes are logically necessary near the top 
and bottom of the PBL, where the organized vertical motions associated with the convective 
circulations must vanish. The time scale for destruction of the convective circulations by lateral 
mixing due to smaller-scale turbulence is T dil . The source of \\f per unit mass per unit time is 
denoted by S ¥ . This could represent, for example, the effects of radiation. 

From (2.1) and (2.2) and the correponding "downdraft" equations, and using (2.3), we 
canrecover the area-averaged conservation equations for mass and y: 

fp = -v.(pv)-pf, 

(2.4) 
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Naturally, all terms involving p, v, and have dropped out of (2.4-5). In (2.5), the turbulent 
fluxes associated with the convective circulations are represented by 


gf ¥ = -coTi' = -[(©. - a>)( Yu ~ Y)° + {o> d - co)( ¥d - y)(l - cr)] 
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is a convective mass flux, which has been defined in such a way that it is less than or equal to zero. 
The total turbulent flux of y is F ¥ = + f ¥ . 
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From the "updraft" mass conservation equation (2.1) and the corresponding "downdraft" 
equation, we can derive an expression for do/di 

do -d d 
~^ + C0 ^ (7 + ^ C0 ' = ^~ V - 

at dp dp 
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By neglecting the local time rate of change and large-scale advection terms, we can simplify (2.8) 
to a more recognizable form: 

d( 0 j dp = fi- v. 


(2.9) 


According to (2.9, the vertical distribution of to. is closely related to the profiles of jj. and v. 

We can also derive the conservation equation for the variance of yrthat is associated with 
the convective circulations: 
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where 
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( 2 . 11 ) 

In (2.10), we have neglected the local time-rate-of-change and large-scale advection terms, 
although we use them in our numerical model. According to (2.12), transport balances the 
combined effects of dissipation and production. Using conventional Reynolds averaging, we find 
corresponding to (2.10) that 
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Here £y is the molecular dissipation rate. Comparison of (2.10) with (2.12) reveals a term-by-term 
correspondence. Note that, in (2.10), the sign of the triple correlation is determined by the value 
of o. Variance is transported upward for o < 1/2, and downward for o > 1/2. 

The damping factor a(l-o) x^, which appears in (2.10), can be simplified by noting that 
small eddies will be most effective at reducing the differences between the ascending and 
descending branches of the plumes when one of the two branches is much narrower than the other, 
and least effective when the two branches have the same width. On this basis, we assume that Tj/ f 
= cj(I-cj) tii is independent of c. We also adopt the definition (T^) 1 = -A co./8p M , where A is 
positive and nondimensional. This is motivated by the idea that the capacity of small eddies to 
dissipate the convective circulations by lateral diffusion increases as co. increases, and as the PBL 
depth decreases, so that A should be less variable than T^. 
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Assuming that the inflow and outflow layers do not overlap, we can eliminate u and v in 
(2.12), to obtain 


where 
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(2.14) 


(2.15) 


The to. term of (2.14) is never negative, so it contributes to X in the same sense as the positive A 
term. It follows that X > 0 unless a increases strongly upward. The 0). and A terms of (2.14) 
both arise from dissipative processes. The Bo/d p term of (2.14) arises from the transport term of 
( 2 . 10 ). 


3. Discussion 

Several conclusions can be drawn from (2.13-14), by considering the case S ¥ = 0, so that 


Zy = 0. 


First, if (1-2 o)/X is small (o near 1/2 and strong dissipation), we get 

m co. d if/ 

= 


(3.1) 


This is a downgradient diffusion formula, with diffusion coefficient -to. / (gX). 

Next, consider a layer that is well-mixed in several conservative variables \|/ 2 , etc. In 
order to satisfy (2. 1 3) simultaneously for all of the variables we must take 


X = 0, 
o = 1/2. 


(3.2) 

(3.3) 


Since, according to (3.3), o is a constant, (2.14) implies that X is the sum of two non-negative 
terms. Then (3.2) implies that each term must be zero, so that 
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c o . = constant, 


(3.4) 


Referring back to (2.10), we find a trivial balance in which gradient production, transport, and 
dissipation are all zero. (When a = 1/2, the variance transport by triple correlations vanishes.) 

Finally, suppose that X is small (weak dissipation), and that a is small. Then (2.14) 
reduces to 


dfv _ 0Kdy_ 

dp ~ g dp 

(3.6) 

This balance is characteristic of the temperature and moisture budgets of tropical cumulus layers, in 
which "compensating subsidence" produces the convective effects on the mean profiles. Arakawa 
(1969) derived (3.6); it has also been discussed by Ooyama , among others. [The "detrainment" 
term included in cumulus parameterization theories is also present in (2.13-14); it arises from the 
gradient of to. in (2.14).] 


For the cumulus regime, the conservation equation for y can often be approximated by a 
balance between vertical advection and convective transport, i.e. 


df ¥ __ co d y 
dp g dp 


Comparing (3.6) and (3.7), we find that 

a = (l - (o, / <o) / 2. 


(3.7) 


(3.8) 


For tropical cumulus layers with large-scale rising motion, to. = to, so that (3.8) is consistent with 
observations that show o « 1. For the case of large-scale sinking motion, (3.8) predicts larger 
values of a, again in qualitative agreement with observations. 

These results show that (2.14-15) are consistent with a variety of observed balances in 
convective layers. They provide a dynamical interpolation between the "compensating subsidence" 
and "mixing length" regimes. 
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